Please type a plus sign (+) inside this box 



PTO/SB/05 (08-00) | 
Approved for use through 9/30/2000. OMB 0651-0032 — — 
Patent and Trademark Office: U.S DEPARTMENT OF COMMERCE I 
Under the Paperwork Reduction Act of 1995, no persons are required to respond to a collection of information unless it displays a valid OMB control number 



UTILITY 
PATENT APPLICATION 
TRANSMITTAL 



Attorney Docket No. | 



7416/78600 



First Inventor or Application identifier \ Bhattacharya, ct al. 

Title | Preprocessing Modules For Quality Enhancement of MBE Codcrsft^ =j 
Express Mail Label No. | EL709889027US ^ 




APPLICATION ELEMENTS 

See MPEP chapter 600 concerning utility patent application 



□ *Fee Transmittal Form (e.g., PTO/SB/17) 
(Submit an original, and a duplicate for fee processing) 

□ Applicant claims small entity status. 
See 37 CFR1.27. 

■■El 



Assistant Commissioner for Patents 3 
ADDRESS TO: Box Patent Application 
Washington, D.C. 20231 

CM 



Specification fj otal Pages 
(preferred arrangement set forth below) 

- Descriptive title of the invention 

- Cross Reference to Related Applications 

- Statement Regarding Fed sponsored R&D 

- Reference to sequence listing, a table, 
or a computer program listing appendix 

- Background of the Invention 

- Brief Summary of the Invention 

- Brief Description of the Drawings (if filed) 

- Detailed Description 

- Claim(s) 

- Abstract of the Disclosure 



35 



7 I I CD-ROM or CD-R in duplicate, large table or ov 

' 1 1 Computer Program (Appendix) H 

8. Nucleotide and/or Amino Acid Sequence Submission 

(if applicable, all necessary) 

a. | | Computer Readable Form (CRF) 

b. Specification Sequence Listing on: 

i. □ CD-ROM or CD-R (2 copies); or 

ii. □ paper 

c. | | Statements verifying identity of above copies 



4. | s/ | Drawing(s) (35 U.S.C. 113) [Total Sheets 



b. 



□ 



Oath or Decla ration [Total Pages 

a. | | Newly executed (original or copy) 

Copy from a prior application (37 CFR 1 -63(d)) 
(for continuation/divisional with Box 17 completed) 

DELETION OF INVENTOR(S) 
Signed statement attached deleting inventor(s) 
named in the prior application, see 37 CFR 
1 .63(d)(2) and 1 .33(b). 

6. | s/ | Application Data Sheet. See 37 CFR 1.76 



□ 



ACCOMPANYING APPLICATION PARTS 



Assignment Papers (cover sheet & document(s)) 

37 CFR 3.73(b) Statement I I Power of Attorri ey 

(when there is an assignee) I 1 

English Translation Document (if applicable) 

Information Disclosure 
Statement (IDS)/PTO-1449 

Preliminary Amendment 



□ Copies of IDS 
Citations 



Return Receipt Postcard (MPEP 503) 
(Should be specifically itemized) 

Certified Copy of Priority Document(s) 
(if foreign priority is claimed) 

unexecuted Declaration and Power of Attorney (4 
pp.); Certificafcof Express Mail (1 p.) 



Other: 



1y the requisite information below and in a preliminary amendment, 



17. If a CONTI NU I N G AP P LI CAT ION, check appropriate box, and suppl 
or in an Application Data Sheet under 37 CFR 1.76: 

Q Continuation Q Divisional Q Continuation-in-part (CIP) of prior application No.: 

Prior application information: Examiner . Group / Art Unit 

For CONTINUATION or DIVISIONAL APPS only : The entire disclosure of the prior application, from which an oath or declaration is supplied under 
Box 5b is considered a part of the disclosure of the accompanying continuation or divisional application and is hereby incorporated by reference. 
The incorporation can only be relied upon when a portion has been inadvertently omitted from the submitted application parts. 



18. CORRESPONDENCE ADDRESS 



fXl Customer Number or Bar Code Label 




or 1 I Correspondence address below 





Walter J. Kawula, Jr. 




24628 








Name 


Welsh & Katz, Ltd. 




PATENT .TRADEMARK OFFICE 








Address 


120 South Riverside Plaza 


22nd Floor 














City 


Chicago 


State 


IL 


Zip Code 


60606 


Country 


USA 


Telephone 


312/655-1500 


Fax 


312/655-1501 



Name (Print/Type) 



Signature 



Walter J. Kawula, Jr. / f\ 

num.. bfouip^ 

This form Isestlmated to take 0.2 houfcg to complete 



Registration No. (Attorney/Agent) 



Date 



39,724 



October 26, 2000 



3 



Burden Hour Statement: This form is estimated to take 0.2 hoi^ to complete. Time will vary depending upon the needs of the individual case Any 
comments on the amount of time you are required to complete this form should be sent to the Chief Information Officer Patent and Trademark Office, 
Washington, DC 20231. DO NOT SEND FEES OR COMPLETED FORMS TO THIS ADDRESS. SEND TO: Assistant Commissioner for Patents, Box 
Patent Application, Washington, DC 20231. 



IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

Applicant: Bhattacharya, et al. 

For: Preprocessing Modules For Quality Enhancement 
of MBE Coders and Decoders For Signals Having 
Transmission Path Characteristics 

Attorney Docket No.: 7416/78600 

Box PATENT APPLICATION 
Commissioner for Patents 
Washington, D.C. 20231 

Sir: 

I hereby certify that the attached: return post card; this Certificate of Express Mail (1 p.); 
Utility Patent Application Transmittal (1 p.); Specification (35 pp.); informal drawings (7 pp.); 
unexecuted Declaration and Power of Attorney (4 pp.); and Application Data Sheet (2 pp.), are 
being deposited with the United States Postal Service via Express Mail No. EL709889027US in an 
envelope addressed to: Box Patent Application, Commissioner for Patents, Washington, D.C. 
20231, on this date. 

October 26, 2000 pJtUL^ 

Rebecca Ruiz 

Welsh & Katz, Ltd. 

120 South Riverside Plaza 

22nd Floor 

Chicago, Illinois 60606 
Phone: 312/655-1500 
Fax: 312/655-1501 




Application Data Sheet 
Inventor Information 

Inventor One Given Name:: 
Family Name:: 
Postal Address Line One:: 
City- 
State or Province:: 
Postal or Zip Code- 
Citizenship Country- 
Inventor One Given Name- 
Family Name:: 
Postal Address Line One- 
City: : 

State or Province:: 
Postal or Zip Code:: 
Citizenship Country- 
Inventor One Given Name:: 
Family Name:: 
Postal Address Line One- 
City:: 

State or Province:: 
Postal or Zip Code:: 
Citizenship Country:: 

Correspondence Information 

Name Line One:: 
Name Line Two:: 
Address Line One:: 
Address Line Two- 
City:: 

State or Province:: 
Postal or Zip Code- 
Telephone Number:: 
Fax:: 



Puranjoy 
Bhattacharya 

203, 8th Main, Mathikere Extension 
Bangalore 

560 054 
INDIA 

Manoj Kumar 
Singhal 

2337, 2nd "A" Cross, 16th "B" Main, Indira Nagar, 

HAL 2nd Stage 

Bangalore 

560 075 
INDIA 



Sangeetha 

1007, 26th "B" Main, 9 th Block, Jayanagar 
Bangalore 



INDIA 



Welsh & Katz, Ltd. 
Walter J. Kawula, Jr. 
22nd Floor 

120 South Riverside Plaza 

Chicago 

IL 

60606 

(312) 655-1500 
(312) 655-1501 



1 



Application Information 



Title Line One:: Preprocessing Modules For Quality Enhancement Of MBE 

Coders And Decoders For Signals Having Transmission 
Path Characteristics 

Total Drawing Sheets:: 7 

Formal Drawings?:: No 

Application Type:: Utility Patent 

Docket Number:: 741 6/78600 



Representative Inforr 

Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 
Registration Number 



One:: 


24,003 


Two:: 


22,839 


Three:: 


28,903 


Four:: 


27,429 


Five:: 


25,060 


Six- 


22,053 


Seven:: 


27,466 


Eight- 


29,434 


Nine:: 


29,054 


Ten:: 


29,381 


Eleven:: 


34,044 


Twelve:: 


27,600 


Thirteen:: 


34,137 


Fourteen:: 


39,724 


Fifteen:: 


37,963 


Sixteen:: 


41,050 


Seventeen:: 


34,177 


Eighteen:: 


47,219 



Domestic Priority Information 



This application claims priority to U.S. Provisional Application 60/161,745, filed 
October 26, 1999 



2 



Utility Patent Application 



RPA 2; 7416/78600 



TITLE OF THE INVENTION 

PREPROCESSING MODULES FOR QUALITY ENHANCEMENT OF MBE CODERS AND 
DECODERS FOR SIGNALS HAVING TRANSMISSION PATH CHARACTERISTICS 

CROSS REFERENCE TO RELATED APPLICATION 

This application claims the benefit of U.S. Provisional Application No. 60/161/745, filed 
October 26, 1999. 

FIELD OF THE INVENTION 
The invention relates to processing a speech signal. In particular, the invention relates to 
enhancing speech signal quality. 

BACKGROUND OF THE INVENTION 

There has been a substantial amount of effort in developing toll-quality speech coders 
that operate below 4kbps. Most of the coders in this bit-range are parametric in nature; One of 
the most prominent among these is the Multiband Excitation (MBE) Coder developed by Griffin 
and Lim, The MBE scheme is derived from mainstream sinusoidal coding (McAulay et al.), 
where voiced speech is reproduced as a weighted sum of sine waves at the harmonics of a pitch 
frequency and unvoiced speech bands are reproduced as bandlimited white noise with 
appropriate amplitudes. The encoding is performed by splitting the input speech into frequency 
bands centered around the harmonics, and recording the respective spectral amplitudes based on 
the outcome of corresponding voicing decisions (assuming the excitation is a sinusoid or 
narrowband noise for the voiced and unvoiced cases, respectively). 

The MBE coding scheme has the potential to produce high quality (in terms of 
intelligibility and naturalness) output speech (Tian et al.) at very low bit rates. The parameters 
used in the MBE coding scheme are also resistant to moderate levels of noise (15 dB wideband 
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white noise). There are, however, some undesirable characteristics of the scheme that severely 
hamper the deployment of MBE-based codecs for the purpose of coding speech produced in 
noisy ambient conditions (above lOdB wideband noise) and/or speech received via transmission 
paths, such as a telephone channel. 

Under transmission path conditions, and in particular, under telephone-channel- 
bandwidth (TCB) conditions, the baseband frequencies are grossly attenuated, as shown in 
Figure 3. This frequently results in the loss of a pitch component and the components first one 
or two harmonics for low-pitched speakers, a phenomenon which greatly hampers pitch 
detection. Pitch detection also becomes increasingly faulty above 15 dB wideband (ambient) 
noise in the input speech signal. However, all parameter estimates in the MBE scheme are, in 
one way or the other, derived via a spectral matching mechanism which in turn crucially depends 
on the harmonic structure created using the pitch parameter. The pitch parameter, therefore is 
the pivotal element in the parameter estimation scheme and, consequently, errors in pitch 
detection frequently lead to the corruption of other parameters. As a result, the MBE codec 
decoded output is prone to several audible distortions such as voice-breaks, screeches, clicks, 
varying levels of hoarseness, and occasional synthetic tonality, for speech having transmission 
path characteristics, such as TCB and/or noisy input speech. 

It has been confirmed, through repeated tests for speech decoded from TCB inputs, that 
voice-breaks observed are frequently associated with pitch region (period) halving, while 
hoarseness is associated with undervoicing. These problems are dominant for low-pitched 
speakers. Tonality, on the other hand, results from overvoicing. 
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One spectral amplitude quantization technique involves intermediate spectral smoothing 
(e.g. if LPC is used, as suggested by Kondoz, a screeching effect is produced for pitch doublings, 
although such occurrences are relatively infrequent). 

The robustness problems discussed above have greatly limited the deployment of MBE 
coders in real-life situations, except for mobile communications, which have significantly lower 
quality demands. In a broad sense, these problems have deterred the achievement of toll-quality 
speech (implying indistinguishable from telephone speech quality) for MBE coders. 

This is unfortunate since MBE coders, which have high compression ratios, may be used 
in a number of applications (primarily storage applications) that are strapped for memory 
resources. The MBE coders provide twice, and in some cases three times the speech storage 
capacity over conventional CELP coders. 

CELP coders imply waveform coding (as opposed to spectral coding in MBE), and 
degrade miserably when operating at rates below 5 kbps. For clean 4kHz bandwidth speech (i.e. 
sampled at 8 kHz, but not subject to the exact telephone-channel frequency response), MBE 
codecs deliver virtually the same output quality, at 2-3 kbps, as higher bit rate (5-6 kbps) CELP 
codecs. However, because of the earlier cited MBE coder problems, the latter continue to be 
preferred for use in voice communication and storage applications that assume noise and 
transmission path characteristics, such as telephone channel bandwidth conditions (CELP codecs 
degrade gracefully under either condition), under normal operating conditions. 

Quality degradation in MBE codecs for noisy and transmission path speech, such as 
telephone-channel speech, has been persistent since its advent. A root cause analysis of the 
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reasons for the distortions induced under the above-mentioned conditions were presented by 
Bhattacharya et ah in 1999, but researchers have been aware of the existence of the problems for 
a long time. 

Researchers, thus far, have attempted to provide robustness to MBE coders by changing 
the basic MBE codec modules. They have essentially suggested alternative methods for the 
robust estimation of pitch and voicing parameters. 

These alternate attempts to compensate for transmission path characteristics, such as 
telephone-channel characteristics, by inverse filtering and to compensate for noise in the input 
signal by spectral subtraction have not been popular mainly because of the associated 
implementation problems. In the former case, designing a stable inverse filter for the telephone 
channel becomes an insurmountable problem when conventional design methods are applied. 
This is because the telephone channel inverse characteristic involves a major gently sloping 
segment accompanied by sharp peaks at either end, and deviation from the expected curve 
becomes audible at virtually all frequencies. In the latter case, the noise compensation process 
breeds a tonal noise called musical noise, which appears at the decoded output as an 
unacceptable distortion. 

Previous solutions to the projected problem have only been marginally effective because 
the basic speech signal is often highly corrupted and because the basic speech signal produces a 
spurious signal with parameter values lying within expected bounds. A common example, in this 
regard, is where a multiple of the pitch frequency becomes the dominant lowest harmonic and 
suppresses the actual fundamental frequency under telephone-channel bandwidth conditions. 
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The amount of parametric corruption varies within wide limits (e.g. depending on the loudness 
and type of noise) farther complicating the robust-estimation process. 

In addition, one should note that there have not been any estimation processes that have 
been 100% reliable even under absolutely clean input speech conditions. The pitch estimation 
accuracy of the invention, when used with the MBE model, decreases gracefully from a 0.2% 
coarse error rate at 30dB ambient (white) noise to a 5% coarse error rate at 10 dB ambient noise. 

Publications relevant to processing signals representing speech include: McAulay et al., 
"Mid-Rate Coding based on a sinusoidal representation of speech", Proc. ICASSP85, pp. 945- 
948, Tampa, Fla., Mar. 26-29, 1985 (discusses the sinusoidal transform speech coder); Griffin, 
"Multi-band Excitation Vocoder", Ph.D. Thesis, M.LT, 1987, (Discusses the Multi-Band 
Excitation (MBE) speech model and an 8000 kbps MBE speech coder); SM. Thesis, M.LT, May 
1988, (discusses a 4800 bps Multi-Band Excitation speech coder); McAulay et al. ? 
"Computationally efficient Sine- Wave Synthesis and its applications to Sinusoidal Transform 
coding", Proc. ICASSP 88, New York, N.Y., pp.370-373, April 1988, (discusses frequency 
domain voiced synthesis); D.W. Griffin, J.S. Lim, "Multi-band Excitation Vocoder," IEEE 
Trans. Acoust., Speech, Signal Processing, vol. 36, pp.1223-1235, August 1988; P. Bhattacharya, 
M. Singhal and Sangeetha, "An analysis of the weaknesses of the MBE coding scheme," IEEE 
international conf. on personal wireless communications, 1999; Tian Wang, Kun Tang, Chonxgi 
Feng "A high quality MBE-LPC-FE Speech coder at 2.4 kbps and 1.2 kbps, Dept. of Electronic 
Engineering, Tsinghua University, Beijing, 100084, P.R. Chinna; Engin Erzin, Arun kumar and 
Allen Gersho "Natural quality variable-rate spectral speech coding below 3.0 kbps, Dept. of 
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Electrical and Computer Eng., University of California, Santa Barbara, CA, 93106 USA; 
INMARSAT M voice codec, Digital voice systems Inc. 1991, version 3.0 August 1991; 
A.M.Kondoz, Digital speech coding for low bit rate communication systems, John Wiley and 
Sons; Telecommunications Industry Association (TIA) "APCO project 25 Vocoder description" 
Version 1.3, July 15, 1993, IS102BABA (discusses 7.2 kbps IMBE speech coder for APCO 
project 25 standard); Telephone transmission quality transmission standards, ITU 
Recommendation p. 48; U.S. Pat. No. 5,081, 681 (discloses MBE random phase synthesis); 
Jayant et al., Digital Coding of Waveforms, Prentice-Hall, 1984, (discussing the speech coding in 
general); U.S Patent No. 4,885,790 (discloses sinusoidal processing method); Makhoul, "A 
mixed-source model for speech compression and synthesis", IEEE (1978), pp. 163-166 ICASS 
P78; Griffin et al. "Signal estimation from modified short-time fourier transform", IEEE 
transactions on Acoustics, speech and signal processing, vol. ASSP-32, No.2, Apr. 1984, pp. 
236-243; Hardwick, "A 4.8 kbps multi-band excitation speech coder", S.M. Thesis, M.I.T., May 
1988; Almeida et al., "Harmonic coding: A low bit rate, good quality speech coding technique," 
IEEE (CH 1746-7/82/000 1684) pp. 1664-1667 (1982); Digital voice systems, Inc. "The DVSI 
IMBE speech compression system," advertising brochure (May 12, 1993); Hardwick et al., "The 
application of the IMBE speech coder to Mobile communications," IEEE (1991), pp. 249-252 
ICASSP 91 May 1991; Portnoff, "Short-time fourier analysis of samples speech", IEEE 
transactions on accoustics, speech and signal processing, vol. ASSP-29, No-3, Jun. 1981, pp. 
324-333; Akaike H. ? "Power spectrum estimation through auto-regressive model fitting," Ann. 
Inst. Statist. Math., Vol. 21, pp. 407-419, 1969; Anderson, T.W., "The statistical analysis of time 
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series" Wiley, 1971; Durbin, J., "The fitting of time-series models," Rev. Inst. Int. Statist., Vol. 
28, pp. 233-243, 1960; Makhoul J., "Linear Prediction: a tutorial review," Proc. IEEE, Vol. 63, 
pp. 561-580, April 1975; Kay S. M., "Modern spectral estimation: theory and application," 
Prentice Hall, 1988; Mohanty M., "Random signals estimation and identification," Van Nostrand 
Reinhold, 1986. The content of the publications listed above are incorporated herein by 
reference. 

BRIEF SUMMARY OF THE INVENTION 

The invention enhances MBE coder performance so that speech having transmission path 
characteristics, such as telephone-channel bandwidth (TCB) and/or noisy speech input, will have 
close to toll-quality speech quality. Pursuant to first and second aspects of the invention, 
separate prefilter and parameter preprocessor modules can be used with an MBE encoder and an 
MBE decoder, respectively. 

Pursuant to a first aspect of the invention, the prefilter module incorporates an inverse 
filter. The effect of the inverse filter compensates for a transmission path transfer function, such 
as a telephone channel transfer function. The frequency domain for a telephone-channel inverse 
filter comprises a smooth middle portion with sudden peakiness at extremities, allowing efficient 
modeling through an all-pole filter. A transfer function of the inverse filter should conform with 
a target characteristic over the entire frequency range (this is in contrast to pass band and stop 
band conventional filters, which have associated gains). The inverse filter can assume the shape 
of an effective all-pole filter and can be of low order, such as, for example, 6 poles. Hence, it is 
computationally efficient. 
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An inverse filter design procedure also ensures that the filter is stable and extremely close 
to desired characteristics. The inverse filter design procedure is general and may be used under 
similar design constraints (i.e. to realize spectra that are peaky or have sudden deep valleys). In 
this case, the inverse characteristic having peaks is used to design an all-pole filter whose 
coefficients are used for an FIR realization of the target spectral characteristic. 

In traditional parametric encoding, it is assumed that corrupted parameters are not subject 
to further improvement. Further, parametric correlation among a series of adjacent frames is 
usually not utilized. Consequently, rectifying encoded parameters for a parametric encoder using 
evolution trajectory information is novel. 

A parameter preprocessor (PP) pursuant to a second aspect of the invention is a module 
that attempts to rectify erroneous estimates of encoded parameters by taking their respective 
evolution trajectories over a succession of frames into account. This module, therefore, 
effectively restores decoded speech quality irrespective of the origin of distortion at the encoder 
input. The parameter preprocessor further assumes simultaneous availability of parameters over 
a sequence of frames, which is common for storage applications. 

The pitch parameter has been identified as the principal indicator of parametric 
corruption at the individual frame level for the MBE coder. Also, since each parameter has been 
found to exhibit characteristic trajectory traits, differing methods have been derived to rectify 
each kind of parameter. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Further objects of the invention, taken together with additional features contributing 
thereto and advantages occurring therefrom, will be apparent from the following description of 
the invention when read in conjunction with the accompanying drawings, wherein: 

Figure 1 depicts a block diagram of a Multiband Excitation encoder that may be used in 
conjunction with the invention; 

Figure 2 depicts a block diagram of a Multiband Excitation decoder that may be used in 

conjunction with the invention; 

Figure 3 depicts the amplitude and frequency characteristics of an IRS filter that models 
spectral characteristics of a telephone channel and that complies with ITU-R (p. 48) 
specifications; 

Figure 4 depicts a block diagram of a prefilter module acting in conjunction with an MBE 
encoder, pursuant to a first aspect of the invention; 

Figure 5 depicts a block diagram of a parameter preprocessor module acting in 
conjunction with an MBE decoder, pursuant to a second aspect of the invention; 

Figure 6 depicts a block diagram of an autoregressive model used to model an inverse 
filter, pursuant to a first aspect of the invention; 

Figure 7 depicts a block diagram of a pitch rectification procedure used in the parameter 
preprocessor module that incorporates principles of a second aspect of the invention; and 

Figure 8 depicts a block diagram of a voicing parameter correction procedure used in the 
parameter preprocessor module that incorporates principles of a second aspect of the invention. 
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DETAILED DESCRIPTION OF THE INVENTION 
While the invention is susceptible to use in various forms and embodiments, there is 
shown in the drawings and will hereinafter be described a specific form and embodiment with 
the understanding that the disclosure is to be considered an exemplification of the invention and 
is not intended to limit the invention to the specific form or embodiment illustrated. 

A block diagram of one MBE encoder that can be used in conjunction with the invention 
is shown in Figure 1 (other encoders not shown may also be used in conjunction with the 
invention). 

The encoder of Figure 1 involves analysis of input speech, parameterization of features 
and quantization of parameters. In the analysis stage of the shown encoder, the input speech is 
passed through block 100 to high-pass filter the signal to improve pitch detection, for situations 
where samples are received through a telephone channel. The output of block 100 is passed to a 
voice activity detection module, block 101. This block performs a first level active speech 
classification, classifying frames as voiced and voiceless. The frames classified voiced by block 
101 are sent to block 102 for coarse pitch estimation. The voiceless frames are passed directly to 
block 105 for spectral amplitude estimation. 

During coarse pitch estimation (block 102) of the encoder shown in Figure 1, a synthetic 
speech spectrum is generated for each pitch period at half sample accuracy, and the synthetic 
spectrum is then compared with the original spectrum. Based on the closeness of the match, an 
appropriate pitch period is selected. The selected coarse pitch is further refined to quarter sample 
accuracy in block 103 by following a procedure similar to the one used in coarse pitch 
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estimation. However, during quarter sample refinement, the deviation is measured only for 
higher frequencies and only for pitch candidates around the coarse pitch. 

In the encoder of Figure 1, based on the pitch estimated in block 103, the current 
spectrum is divided into bands and a voiced/unvoiced decision is made for each band of 
harmonics in block 104 (a single band comprises three harmonics). For each band, a spectrum is 
synthesized first assuming all the harmonics in the band are voiced and then assuming all the 
harmonics in the band are unvoiced. An error for each synthesized spectra is obtained by 
comparing the respective synthesized spectrum with the original spectrum over each band. If the 
voiced error is less than the unvoiced error, the band is marked voiced, otherwise it is marked 
unvoiced. 

In the encoder of Figure 1, a Voicing Parameter (VP) is introduced to reduce the number 
of bits required to transmit the voicing decisions found in block 104. The VP denotes the band 
threshold, under which all bands are declared unvoiced and above which all bands are marked 
voiced. Instead of a set of decisions, a single VP is calculated in block 107. 

Speech spectral amplitudes are estimated by generating a synthetic speech spectrum and 
comparing it with the original spectrum over a frame. The synthetic speech spectrum of a frame 
is generated so that distortion between the synthetic spectrum and the original spectrum is 
minimized in a sub-optimal manner in block 105. 

Spectral magnitudes are computed differently for voiced and unvoiced harmonics. 
Unvoiced harmonics are represented by the root mean square value of speech in each unvoiced 
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harmonic frequency region. Voiced harmonics, on the other hand, are represented by synthetic 
harmonic amplitudes, which characterize the original spectral envelope for voiced speech. 

The spectral envelope contains magnitudes of each harmonic present in the frame. 
Encoding these amplitudes require a large number of bits. Because the number of harmonics 
depends on the fundamental frequency, the number of spectral amplitudes varies from frame to 
frame. Consequently, in the encoder of Figure 1, the spectrum is quantized assuming it is 
independent of fundamental frequency, and modeled using a linear prediction technique in 
blocks 106 and 108. This helps reduce the number of bits required to represent the spectral 
amplitudes. LP coefficients are then mapped to corresponding Line Spectral Pairs (LSP) in 
block 109, which are then quantized using multi-stage vector quantization. The residual of each 
stage is quantized in a subsequent stage in block 110. 

A block diagram of an MBE decoder that may be used with the invention is illustrated in 
Figure 2 (other decoders not shown may also be used in conjunction with the invention). 

Parameters from the encoder are first decoded in block 200. A synthetic speech spectrum 
is then reconstructed using decoded parameters, including fundamental frequency values, 
spectral envelope information and voiced/unvoiced characteristics of the harmonics. Speech 
synthesis is performed differently for voiced and unvoiced components and consequently 
depends on the voiced/unvoiced decision of each band. Voiced portions are synthesized in the 
time domain whereas unvoiced portions are synthesized in the frequency domain. 

In the decoder of Figure 2, the spectral shape vector (SSV) is determined by performing a 
LSF to LPC conversion in block 201. Then using the LPC gain and LPC values computed 
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during the LSF to LPC conversion (block 201), a SSV is computed in block 202. The SSV is 
spectrally enhanced in block 203 and inputted into block 204. The pitch and VP from the 
decoded stream are also inputted into block 204. In block 204, based on the voiced/unvoiced 
decision, a voiced or unvoiced synthesis is carried out in blocks 206 or 205, respectively. 

An unvoiced component of speech is generated from harmonics that are declared 
unvoiced. Spectral magnitudes of these harmonics are each allotted a random phase generated 
by using a random phase generator to form a modified noise spectrum. The inverse transform of 
the modified spectrum corresponds to an unvoiced part of the speech. 

Voiced speech represented by individual harmonics in the frequency domain is 
synthesized using sinusoidal waves. The sinusoidal waves are defined by their amplitude, 
frequency and phase, which were assigned to each harmonic in the voiced region. 

The phase information of the harmonics is not conveyed to the decoder. Therefore, in the 
decoder of Figure 2, at transitions from an unvoiced to a voiced frame, a fixed set of initial 
phases having a set pattern is used. Continuity of the phases is then maintained over the frames. 
In order to prevent discontinuities at edges of the frame, due to variations in the parameters of 
adjacent frames, both the current and previous frame's parameters are considered. This ensures 
smooth transitions at boundaries. The two components are then finally combined to produce a 
complete speech signal by conversion into PCM samples in block 207. 

Pursuant to first and second aspects of the invention, separate prefilter and parameter 
preprocessor modules are used with an encoder, such as, for example, the MBE encoder depicted 
in Figure 1, and a decoder, such as, for example, the MBE decoder depicted in Figure 2, 
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respectively. As a result, on one hand, the modules do not preclude further developments to the 
MBE coder structure, and on the other, the modules may be strapped onto various 
implementations of the MBE coder, including hardware implementations. 

Two modules may be used, one for preprocessing the input signal before it enters the 
encoding process (Figure 1), and the other for preprocessing encoded parameters before they are 
processed by the decoder (Figure 2). These modules will be referred to as the prefilter and 
parameter preprocessor (PP) modules respectively. Either of these can operate in isolation of the 
actual MBE codec modules. Consequently, an improvement to the basic MBE models 
necessarily accrue to the augmented configuration. 

Figure 4 shows a block-level signal flow through a prefilter module and MBE encoder, 
pursuant to a first aspect of the invention. The input signal (Block-4.1) is processed by the 
prefilter module (Block-4.2) to produce a transmission path compensated signal, and in 
particular, a telephone-channel-compensated signal. The compensated signal is encoded by a 
MBE encoder (Block-4.3) to produce the encoded parameter stream (Block-4.4). 

Figure 5 shows a block-level signal flow through a parameter pre-processor and MBE 
decoder, pursuant to a second aspect of the invention. The input encoded parameter stream 
(Block-5.1) is processed by the parameter preprocessor module (Block-5.2) to produce an error- 
corrected parameter stream, which is subsequently decoded by a MBE decoder (Block-5.3) to 
produce output speech (Block-5.4). 

Pursuant to a first aspect of the invention, the prefilter module used in conjunction with 
an MBE encoder incorporates an inverse filter. The inverse filter can be designed to preprocess 
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input speech that has transmission path characteristics, such as TCB speech, by restoring the 60- 
200 Hz band eliminated during transmission through telephone channels. One type of inverse 
filter pursuant to a first aspect of the invention comprises an all-pole filter that can be strapped 
on to the input stage of a MBE speech encoder. 

The inverse filter may be characterized as having an inverse amplitude characteristic of 
the amplitude characteristics of an IRS filter (details in ITU-R P. 48, shown in Figure 3) or other 
filters that approximate frequency-amplitude characteristics of a transmission channel. The IRS 
filter approximates frequency-amplitude characteristics of a telephone channel. In other words, 
if the IRS filter has a frequency response H{co) and the inverse filter has a frequency 
response G(co) , then G(w) should be characterized by the following relation: 

\H(a>)\ 

The desired inverse characteristic of the filter has extremely sharp transitions around 200 
Hz and 3300 Hz, further, the intermediate region has a variable slope. As a result, FIR or IIR 
filters designed by available procedures are lacking. 

It should be noted that an all-pole filter is well suited in the context of an inverse filter 
because of an all pole filter's capability to fit peaky spectral characteristics, and therefore an 
inverse filter solution within this restricted class of IIR filters is beneficial. An inverse filter, 
illustrated below, is one example of such an all-pole filter. One method to design the illustrated 
inverse filter using spectral estimation theory is described below. 
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In this disclosure, the IRS filter is described by the function h(t) in the time domain and 
the illustrated inverse filter is described by the function g(t) where H(co) is the Fourier transform 
of h(t) and G(co) is the Fourier transform of g(t). The objective is to design the illustrated inverse 
filter so that 

G(o>)« (2) 

\H(co)\ 

One method of meeting the objective is to represent a random signal with a power 
spectral density (PSD) equal to | G(a>) | 2 using an auto-regressive (AR) model. An AR model 
that comprises an all-pole system excited by white noise e(n), as shown in Block 6.1 of Figure 6, 
can be used to represent | G(co) | 2 . 

The output sequence of Figure 6 [g(n)] may be characterized by the equation: 

p 

g{n) = -Y j a k S{n-k) + e{n) (3) 

where e(n) is an additive white Gaussian noise sequence. The white noise e(n) has a unit power 
spectral density by definition and the PSD of the random signal being modeled is equal to the 
square of the magnitude response of the all-pole filter. 

Substituting phase information of the inverse filter by a random sequence g(n) allows the 
above described transformation. Note that this transformation is possible because a phase 
characteristic restriction of an inverse filter has not been imposed. In addition, note that the 
assumed random phase is never explicitly specified or used in the design process. 
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iG(G>)r= 



The power spectral density of G(co) may be characterized by the equation: 
1 



(4) 



|l + 2>* 



-Jojk 1 2 



The parameters of an AR model (ak) can be obtained from the auto-correlation function 
(ACF) of the random signal by setting up Yule- Walker equations as follows: 



R(0) R(-Y) - R(-p) 
R(l) R(0) 

R(p) R(p-l) - *(0). 



" 1 " 
















0 






0 


"p. 




_ 0 _ 



(5) 



where R(i) = R(-i), i = l,... 9 p , are the respective ACFs at various lags, and a 2 p is the 
"minimum mean-squared prediction error" for the AR model, which is also equal to the variance 
of the assumed input white noise sequence. 

The ACF R(m) of the virtual random signal g(n) employed in the above equations can 
be efficiently estimated as the inverse Fourier transform of its PSD (Wiener-Khintchine 
Theorem), which, under the given circumstances is equal to the square of the inverse magnitude 
characteristic. This is characterized by the following equation: 

jlTtkm 



R(m) = Y ^— r-e 



(6) 



k=0 | 



The Yule-Walker equations can be solved using a variety of methods, including the 
Levinson-Durbin algorithm which exploits the Toeplitz structure of the leftmost matrix in 
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equation 5. The coefficients ( a x , . . . , a p ) of equation (5) are solved for and used to determine the 

illustrated inverse filter, which is one example of a suitable all-pole filter. 

The illustrated inverse filter may be designed using several methods, the following steps 
describe one method to design the illustrated inverse filter. 

1. Assume the IRS filter is specified as a sequence h(n), n=0,l,...,N-l. 

2. Obtain a new sequence h x (n) by padding zeroes to make the sequence length equal to the 
nearest power of 2, say M. 

3. Obtain H x (k) 9 k=0,l . .,M-1 as the Fast Fourier Transitions of the sequence h x («) , 
n=0,l,...,M-l. 

4. Obtain P(k) = — - , k=0,l,...,M-l. 

\H(k)\ 2 

5. Produce R(m) , m=0,l,. . JM-l, by taking the IFFT of the sequence P(k), k=0,l 5 . . .,M-1 . 

6. Set up the Yule- Walker equations, using R(m) computed in step 5, as per equation 5. 

7. Solve Yule-Walker equations produced by assuming a "q" order AR model through the 
Levinson-Durbin method to obtain the required all-pole filter coefficients. 

Those of ordinary skill in the art will note that Step 7 merely requires a solution of the 
Yule- Walker equations, and is amenable to methods other than the Levinson-Durbin method. 

Those of ordinary skill in the art will also note that there are several methods to meet the 
objective of designing the inverse filter. A second method to meet the objective involves 
modeling | H{co) \ 2 using a Moving Average (MA) model. The MA model parameters are found 
by solving a set of equations set up using the Inverse Fast Fourier Transform of | H(a>) \ 2 . These 
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MA model parameters correspond to the numerator polynomial of the direct system, hence they 
also correspond to the denominator polynomial of the desired inverse characteristic, and hence 
are the coefficients of the target IIR filter. The MA parameter estimation problem (frequently 
handled, as mentioned by Kay, through conversion of the MA process into an equivalent AR 
process), lacks a direct computational solution, reducing the viability of the second method. 

In an experiment performed using 15000 frames of telephone-quality test data, the above 
construct was found to eliminate approximately 80% of the audible artifacts for the MBE codec. 
The invention has been rigorously tested in lab, using simulated, as well as actual telephone 
speech data. 

In spite of the efficacy of the tested inverse filter, some audible artifacts may persist. 
Most of these result from erroneous pitch parameter detection as a multiple or sub-multiple of 
the true pitch parameter value. This is caused by, in certain situations, pitch component 
attenuation. For example, when pitch components attenuation occurs other harmonics or sub- 
harmonics may dominate, and these harmonic or sub-harmonics may ultimately be preferred 
during the matching procedure over the true value. These audible distortions can be eliminated 
prior to decoding, for applications (primarily storage applications) by parameter preprocessing a 
parameter stream from the encoder over a succession of frames. 

As discussed earlier, the corruption of various parameter estimates for the MBE model is 
rooted in gross errors in pitch estimation. Pitch parameter corruption, therefore, is used as the 
primary indicator of parameter corruption over individual frames. The first major step in 
parameter preprocessing, therefore, is detecting pitch parameter corruption. 
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The theory behind parameter error detection as well as parameter error correction is 
based on the gradual variation of most parameters (excluding voicing boundaries) over a 
sequence of frames. Consequently, the value of a parameter over a frame may be predicted from 
neighboring parameter values. Pursuant to a second aspect of the invention, the theory of 
gradual variation of parameters over successive frames is utilized to preprocess signal data. 

One example of using the gradual variation involves parameter preprocessing. Parameter 
preprocessing involves correcting gross pitch errors (primarily doubling and halving errors) 
using trajectory information and updating other coded parameters accordingly. For example, one 
method of parameter preprocessing that involves three stages is described below. A first step 
involves pitch rectification, a second step involves updating spectral amplitudes and a third step 
involves updating voicing parameters. 

The first step of parameter preprocessing in the described method involves pitch 
rectification. During real-time operation of the encoder, spectral matching schemes concentrate 
on information contained within the same frame, with minor augmentation using interframe 
dependencies during tracking. In close temporal proximity to the storage phase (i.e. preceding or 
succeeding storage), however, the entire pitch trajectories may be available, and these may be 
processed using continuity constraints because the pitch parameter changes smoothly over 
contiguous (voiced) stretches. Two important tools in this regard are: (1) a linear low-pass filter 
for smoothing, and (2) a median filter. The latter family of filters is efficient for removing 
sudden departures from the trajectory, while the former smoothes the trajectories. In the 
described preprocessing method, a long-order median filter may be followed by a smaller-order 
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smoothing filter to remove a large number of pitch halving and doublings, especially ones that 
occur in smaller chunks (2-3) frames. The filters may be turned off at voiced-region boundaries 
marked by three or more successively occurring unvoiced frames (a voicing parameter maybe 
used to derive voicing information). 

In the described method, the pitch correction procedure involves predicting pitch value 
using the linear and median filters described above. The closest multiple or sub-multiple of the 
actual reported value of P (e.g. 2P, 3P, P/2, P/3 etc.) to the pitch value of the linear and median 
filtered pitch trajectory is selected as the corrected pitch valve. In actual implementations, these 
four derived pitch values are used for comparison, since the possibility of higher multiples and 
sub-multiples occurring is minimal. Those skilled in the art will recognize, however, that any 
number of sub-multiples and/or multiples may be used while selecting a corrected pitch value. 

Figure 7 shows a schematic diagram of the first step (pitch rectification) of the described 
method. The sequence of pitch values is first median filtered (Block-7.1), and then linearly 
smoothed (Block-7.2). The resulting value is then matched to various multiples and sub- 
multiples of the actual reported pitch value (Block-7.3). The closest multiple or sub-multiple 
match is declared as the corrected pitch value. 

As mentioned earlier, mere correction of the pitch value does not automatically rectify 
other respective artifacts because, apart from leading to the proliferation of fine parametric 
errors, the entire banding structure is changed (e.g. when a pitch-period halving occurs, there are 
half as many spectral coefficients recorded). An updating procedure for other parameters, 
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operating over frames with pitch errors, requires band-structure restoration as well as correction 
of minor errors through trend information. 

The second step of parameter preprocessing in the described method involves updating 
spectral amplitudes. In the second step, all pitch errors (gross ones) are classified into halvings, 
doublings, triplings etc. If the pitch frequency originally detected was half the corrected value, 
there will be twice as many harmonics. If a spectrum is reconstructed by deleting odd 
harmonics, the original spectrum will be restored. 

If, on the other hand, the pitch frequency detected originally was twice the corrected 
value, the alternate harmonics have not been computed (i.e. spectral amplitudes). These can, 
however, be partially reconstructed, assuming smoothness of the gross spectrum, by log-linear 
interpolation between alternate harmonics over the same frame. 

Similar schemes of spectral amplitude restoration can be employed for other harmonics 
and sub-harmonics of incorrectly detected pitch frequency. Procedures to modify spectra 
relating to pitch frequencies that were 1/2, 1/3, 2 times, or 3 times the corrected pitch value are 
listed below. Those skilled in the art will recognize that similar procedures may be used to 
modify other spectra. 

For example, if the pitch frequency originally detected was one-half of the corrected pitch 
value, only 2kth harmonics (i.e. the second, fourth, sixth, etc. harmonic) should be retained. If 
the pitch frequency originally detected was one-third of the corrected pitch value, only 3kth 
harmonics (i.e. the third, sixth, ninth, etc. harmonic) should be retained. If the pitch frequency 
originally detected was twice the corrected pitch value one harmonic should be inserted at the 
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(k + !/ 2 )th harmonic position between successive harmonics (i.e., insert a Yik harmonic between 
the 0 and 1 st harmonics, insert a VA k harmonic between the 1 st and 2 nd harmonics, etc). The 
amplitude of the inserted (k + '/ 2 )th harmonic can be characterized by the equation: 

A(k+l/2) = y}A(k)-A(k + V) (7) 
If the pitch frequency originally was three times the corrected pitch value, two harmonics should 
be inserted at (k + l/3)th and (k + 2/3)th positions between successive harmonics (i.e., insert a 
l/3k and 2/3k harmonic between the 0 and 1 st harmonics, insert a 1 l/3k harmonic and 1 2/3k 
harmonic between 1 st and 2 nd harmonics, etc). The amplitudes of the inserted (k + l/3)th and (k 
+ 2/3)th harmonics can be characterized by the equations: 



A(k + l/3) = \lA\k)A(k+l) ( 8 ) 



A(k + 2/3) = \lA(k)A 2 (k + V) 

The third step of parameter preprocessing in the described method involves updating 
voicing parameters. Trajectories of voicing are characterizable during a single voiced-to- 
unvoiced transition, and a Voicing Parameter (VP) is assumed for the spectrum of each frame of 
voiced speech. When the pitch is detected inaccurately, the VP, which is estimated using the 
same spectral matching scheme as the pitch parameter is estimated with, usually plunges 
abruptly to a low value. This, apart from certain extreme cases, does not usually cause the entire 
frame to be detected as unvoiced, therefore preventing circularity in the error correction 
procedure (note that the pitch correction is based on a frame voicing decision derived from the 
VP). 
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Pursuant to the third step of the described method, the VP can be partially restored by 
obtaining an estimate through smoothing a VP trajectory over a small sequence of frames 
centered around the erroneously coded frame (characterized by a detected gross pitch error) 
using median and linear filtering. The filtered value can then be recorded as the corrected VP. 

Figure 8 shows a schematic diagram of the third step (voicing parameter updation) of the 
described method. The input VP sequence is first median filtered (Block-8.1) and subsequently 
linear filtered (Block 8.2) to generate the output VP sequence. 

The described inverse filter and parameter preprocessor were tested using a 15,000 frame 
test sequence. The test showed that the described inverse filter and parameter preprocessor 
minimized observable errors of the 15,000 test frame sequence to levels close to non-TCB (clean 
input speech) levels. In addition, at the expense of a short initial delay, the test showed that the 
described inverse filter and parameter preprocessor can be applied to real time encode-decode 
applications. 

The described error correction procedures operate under the assumption that parameter 
trajectories obtained over frame sequences are reflective of the principal variational trends, and 
that they do not explicitly depend upon the mechanism causing the errors. Therefore, the 
methods for parameter correction through preprocessing are equally applicable to parameter 
degradation in TCB conditions and high levels of input ambient noise. 

From the foregoing it will be observed that numerous modifications and variations can be 
effectuated without departing from the true spirit and scope of the invention. It is to be 
understood that no limitation with respect to the specific use illustrated is intended or should be 
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inferred. The disclosure is intended to cover by the appended claims all such modifications as 
fall within the scope of the claims. 
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CLAIMS 

What is claimed is: 

1 . A method of signal processing signals having transmission path characteristics, 
comprising the steps of: 

inverse filtering an input signal having transmission path characteristics before 
processing the input signal wherein the transmission path characteristics of the input signal are 
reduced; and 

processing the input signal. 

2. The method of claim 1 , wherein an inverse filter is used to filter the input signal 
and an encoder is used to process the input signal, the inverse filter being in communication with 
the encoder. 

3. The method of claim 2 7 wherein the encoder is a Multi-Band Excitation (MBE) 
encoder and the inverse filter produces a frequency response having a smooth middle portion 
with peakiness at extremities of the frequency response, 

4. The method of claim 2, wherein the inverse filter comprises an all pole filter. 

5. The method of claim 4, wherein the inverse filter is a low order filter including 
about six poles. 

6. The method of claim 2 9 wherein the inverse filter has an inverse amplitude 
response of a filter described by h(t) 9 the filter approximating noisy ambient conditions including 
telephone-channel-bandwidth conditions and the pre-filter response being characterized by: 

G(g>) ~ — wherein H(co) is the frequency response oih{f) and G(a>) is the inverse filter 

\H(co)\ 

frequency response. 
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7. The method of claim 6, wherein a random signal having a power density 
characterized by | G(co) \ 2 is used to design the inverse filter, wherein the processing step 
comprises the sub-steps of: 

parameterizing the input signal; and 

encoding the input signal; and the processing the signal method further comprises 

the steps of: 

preprocessing the encoded signal; and 

decoding the preprocessed encoded signal, wherein a parameter preprocessor is 
used to preprocess the encoded signal and a decoder is used to decode the preprocessed encoded 
signal, the encoder being in communication with the parameter preprocessor and the parameter 
preprocessor being in communication with the decoder. 

8. A method for preprocessing a signal having transmission path characteristics, 
comprising the steps of: 

obtaining at least one sequence, wherein one of the at least one obtained sequence 
is a first sequence [h(n)] wherein n = 0,1, . . . N-l, and N-l is a length value of the first sequence; 

taking a Fast Fourier Transform (FFT) of an obtained sequence to determine H(k); 

obtaining P(k) by using H(k), wherein P{k) is characterized by: P{k) ~ — - , 

\H{k)\ 

k=0,l,.. .,M-1, wherein M is a length value; 

taking an inverse Fast Fourier Transform (IFFT) of P(k) to obtain R(m) 9 wherein 
m = 0, 1, ...M-l; 

preparing Yule- Walker equations using the obtained R(m) values; 

solving the Yule-Walker equations to obtain coefficients; 

using the obtained coefficients to design an inverse filter; and 

preprocessing the signal having transmission path characteristics with the inverse 

filter. 
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9. The method of claim 8 wherein the using the obtained coefficients to design the 
inverse filter comprises the sub-steps of: 

using the obtained coefficients to determine G{co% wherein G(co) is a frequency 
response of the inverse filter, and wherein G{co) * ) , H{oi) being the frequency response 

of h{t\ h(t) being a time domain description of a filter that approximates transmission path 
characteristics including telephone-channel-bandwidth conditions, and h(n) being a sequence 
representing the approximating filter; 

using GO) to determine g(t), wherein g(t) is the time domain description of the 
inverse filter; and 

using g{i) to design the inverse filter. 

10. The method of claim 9 wherein \G(a>)\ 2 is characterized by the equation 

\G(co)\ 2 = -J 

ii+5>**-'* i 2 

wherein a k are the obtained coefficients. 

1 1 . The method of claim 8 wherein at least two sequences are obtained in the 
obtaining at least one sequence step, wherein a second sequence [hi(n)] that modifies the first 
obtained sequence [/?(«)] is one of the at least two obtained sequences, the second sequence 
having a length M and the M length value being equal to a closest power of 2 after the N-l 
length value; and 

wherein the FFT is taken on the second obtained sequence [hi(n)] to determine 

H(k). 

12. The method of claim 8 wherein the Yule- Walker equations are characterized by: 
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wherein <j 2 p is a minimum mean-squared error of an auto recursive model, and , . . . , a p are the 
coefficients to be solved for. 

13. The method of claim 8 wherein the inverse Fast Fourier Transform (IFFT) of P{k) 
is characterized by: 

K -i jlnkm 

IFFT(P(£)) = R(m) - J 7—7—7 * e - . 

14. The method of claim 8 wherein the Yule Walker equations are solved by 
assuming a 'q' order Auto Recursive model and using a Levinson-Durbin method. 



15. A method of preprocessing a signal having transmission path characteristics, 
comprising the steps of: 

obtaining a frequency response [H(a>)] of a filter that approximates noisy ambient 
conditions including telephone-channel-bandwidth conditions; 

modeling \H(ai)\ 2 using a moving average model comprising the sub-steps of: 

taking the inverse Fast Fourier Transform (IFFT) of \H(co)\ 2 to formulate a 

set of equations; 

solving the set of equations to obtain moving average model parameters; 
using the moving average model parameters to design an inverse filter; and 
preprocessing the signal having transmission path characteristics with the inverse 

filter. 
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16. The method of claim 15 wherein the using the moving average model average 
model parameters to design the inverse filter comprises the sub-steps of: 
applying the parameters to the equation: 

|G(*;)| 2 = - 



! 1 + |>* 



- jcck 1 2 



\H(a>)\ 



wherein G(a>) is the frequency response of the pre-filter and a k are the model parameters; and 
using G{co) to design the inverse filter. 

17. A method of processing received encoded data, comprising the steps of: 
preprocessoring the received encoded data before decoding the data; and 
decoding the data. 

18. The method of claim 1 7 wherein a parameter preprocessor performs the 
preprocessoring step and a decoder performs the decoding step, the parameter preprocessor being 
in communication with the decoder. 



19. The method of claim 18 wherein the received encoded data is received from an 
encoder, the encoder being in communication with the parameter preprocessor. 

20. The method of claim 17, wherein the preprocessoring the received encoded data 
step includes the sub-steps of: 

obtaining signal data from the received encoded data, wherein the obtained data 
includes pitch parameter data for a trajectory of successive frames of the signal; 

removing at least one pitch parameter departure from the trajectory of successive 

frames; 

smoothing the trajectory; 
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calculating at least one multiple corresponding to an obtained pitch parameter of a 
frame having a pitch parameter departure and at least one sub-multiple corresponding to the 
obtained pitch parameter; 

comparing a pitch parameter from the removed and smoothened trajectory that 
corresponds to the obtained pitch parameter with the at least one corresponding multiple and the 
at least one corresponding sub-multiple; and 

replacing the obtained pitch parameter with a new pitch parameter based on the 
comparison, the new pitch parameter being selected from the at least one corresponding multiple 
and the at least one corresponding sub-multiple. 

21. The method of claim 20, wherein two multiples are calculated for each obtained 
pitch parameter of a frame having a pitch departure, two sub-multiples are calculated for each 
obtained pitch parameter of a frame having a pitch departure, the comparing step is performed 
for each obtained pitch parameter of a frame having a pitch departure, and the replacing step is 
performed for each obtained pitch parameter of a frame having a pitch departure. 

22. The method of claim 20, wherein a medium filter is used to remove departures 
and a small order linear filter is used to smoothen. 

23. The method of claim 22, wherein the obtained data further includes unvoiced 
frame and voiced frame information, and wherein the medium filter and the small order linear 
filter are turned off when three successive unvoiced frames are detected. 

24. The method of claim 20, wherein the obtained data further includes spectral 
amplitude information; and wherein the preprocessoring the received encoded data step further 
includes the sub-steps of: 

adjusting a number of harmonics for a spectrum of a frame having a new 

pitch parameter. 
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25. The method of claim 24, wherein the adjusting a number of harmonics step 
includes the sub-steps of: 

removing each (2k-l)th harmonic of the spectrum if the new pitch parameter is 
one-half the value of the obtained pitch parameter; 

removing each (3k-l)th harmonic and each (3k-2)th harmonic of the spectrum if 
the new pitch parameter is one-third the value of the obtained pitch parameter; 

inserting one harmonic at each (k + 1/2) location of the spectrum if the new pitch 
parameter is twice the value of the obtained pitch parameter, each inserted (k+l/2)th harmonic 
having an amplitude characterized by the equation A (k+ 1/2) = yjA(k) * A(k + 1) ; and 

inserting one harmonic at each (k+1/3) and one harmonic at each (k+2/3) location 
of the spectrum if the new pitch parameter is three times the value of the obtained pitch 
parameter, each inserted (k+l/3)th harmonic having an amplitude characterized by the equation 

A(k + 1/3) = ljA 2 (k)A(k+l) and each inserted (k+2/3)th harmonic having an amplitude 
characterized by the equation A(k + 2 / 3) = %jA{k)A 2 (k + 1) . 

26. The method of claim 24, wherein the obtained data further includes voice 
parameter information; and 

wherein the preprocessoring the received data step further includes the sub-steps 

of: 

medium filtering a voice parameter trajectory, the voice parameter 
trajectory including voice parameter information of the frame having a new pitch parameter, 
voice parameter information of frames preceding the frame having a new pitch parameter, and 
voice parameter information of frames succeeding the frame having a new pitch parameter; 

linear filtering the voice parameter trajectory; 

using the medium and linear filtered voice parameter trajectory to obtain a 
new voice parameter trajectory. 
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27. A speech system comprising: 

an inverse filtering means for inverse filtering signal data having transmission 
path characteristics; and 

an encoder, the encoder including parameterizing means for parameterizing the 
signal data and encoding means for encoding the signal data, the encoder being in 
communication with the inverse filtering means. 

28. The speech system of claim 27 further comprising: 

a parameter preprocessor, the parameter preprocessor including receiving means 
for receiving the encoded signal data and preprocessoring means for preprocessoring the 
received encoded signal data, the parameter preprocessor being in communication with the 
encoder; and 

a decoder, the decoder including decoding means for decoding the preprocessed 
signal data and synthesizing means for synthesizing the preprocessed signal data into a speech 
signal, the decoder being in communication with the parameter preprocessor. 

29. The speech system of claim 28 wherein the encoder further includes storage 
means to store encoded data and transmission means to transmit encoded data. 

30. A speech system comprising: 

an inverse filtering means for inverse filtering signal data having transmission 
path characteristics; 

an encoder, the encoder including parameterizing means for parameterizing the 
signal data and encoding means for encoding the signal data, the encoder being in 
communication with the inverse filtering means; 

a parameter preprocessor, the parameter preprocessor including receiving means 
for receiving the encoded signal data and preprocessoring means for preprocessoring the 
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received encoded signal data, the parameter preprocessor being in communication with the 
encoder; 

a decoder, the decoder including decoding means for decoding the preprocessed 
signal data and synthesizing means for synthesizing the preprocessed signal data into a speech 
signal, the decoder being in communication with the parameter preprocessor. 
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ABSTRACT OF THE DISCLOSURE 
Pursuant to one aspect of the invention, a prefilter module that incorporates an inverse 
filter is used in conjunction with an encoder. The inverse filter has an inverse frequency 
response of a frequency response of a filter that simulates speech having transmission path 
characteristics, such as telephone-channel bandwidth speech, and/or noisy speech. The inverse 
filter is used to compensate transmission path characteristics of an input signal. The inverse 
filter can be designed using several methods, such as, for example, an autoregressive model or a 
moving average model. Pursuant to a second aspect of the invention, a parameter preprocessor is 
used in conjunction with a decoder. The parameter preprocessor performs pitch rectification 
through use of a medium and linear filter, and updates spectral amplitudes and voicing parameter 
depending on the pitch rectification. The inverse filter and parameter preprocessor, used in 
conjunction with an encoder and decoder, respectively, improve signal processing and parameter 
estimation. 
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